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A  STUDY  OF  TECHNOLOGY  FOR  DETECTION  OF  NONLUMINOUS 
ARTIFICIAL  SATELLITES  IN  DAYLIGHT 

(Romanized  title:  "Baitian  Feizifaguang  Tianti  Mubiao  Tance  Jishu  De  Yanjiu ") 

by  He  Cheng  and  Zhu  Qixiang  (Academia  Sinica,  Photoelectric 
Technology  Institute,  Chengdu  610209) 

Abstract:  In  this  paper,  the  authors  discuss  the  detection  limits  of  systems  for  detecting 
nonluminous  artificial  satellites  in  daylight,  as  well  as  the  relationship  of  target  contrast  and 
signal-to-noise  ratio  to  artificial  satellite  characteristics,  sky  background  characteristics,  optical 
system  parameters,  and  CCD  [charge-coupled  device]  parameters.  The  authors  discuss  in  depth 
the  effects  of  spectral  filtering  technology  on  detection  limits,  contrast,  and  signal-to-noise  ratio. 
Results  of  numerical  simulations  demonstrate  that  spectral  filtering  technology  can  increase  the 
target  satellite’s  optical  contrast  by  4  to  8  times,  raise  the  target  satellite’s  signal-to-noise  ratio 
by  1.5  to  2.5  times,  and  enhance  the  system’s  limit  detection  brightness  by  0.4  to  0.7  visual 
magnitude. 

Foreword 

Detection  of  signals  from  low-contrast,  weak  targets  that  reflect  sunlight  in  daytime  when 
background  signals  are  strong  is  a  key  technology  for  optical  and  photoelectric  detection, 
tracking,  and  measurement  of  nonluminous  artificial  satellites.  During  the  daytime,  when  the  sky 
is  bright,  background  signals  are  strong  and  artificial  satellites  shine  by  reflecting  sunlight.  The 
spectral  characteristics  of  artificial  satellites  are  similar  to  those  of  their  backgrounds.  Thus,  it 
is  necessary  first  to  solve  the  technical  problems  of  extracting  the  signals  of  artificial  satellites 
from  intense  background  signals  when  carrying  out  detection  during  the  daytime.  The  methods 
that  are  often  used  today  to  solve  this  problem  include: 

(1)  Adopting  refrigeration  technology  to  reduce  the  device’s  thermal  noise; 
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(2)  Using  a  low-noise,  wide-band  pre-amplifier  to  amplify  the  photoelectric  detector’s 
output  signal.  This  method  can  only  prevent  noise  from  expanding,  and  has  no  suppressing  effect 
on  the  noise  already  produced  by  the  detector; 

(3)  Employing  relevant  techniques  to  suppress  noise  and  increase  the  output  signal-to- 
noise  ratio.  Using  correlated  dual-sampling  circuits  lowers  the  charge-coupled  device’s  reset 
noise,  but  cannot  lower  its  shot  noise; 

(4)  Adopting  self-adaptive  background  suppression  methods  to  suppress  sky  background 
signals.  These  methods  cannot  reduce  the  shot  noise  created  by  background  photon  fluctuation, 
either. 


In  this  paper,  we  advance  an  optical  method  of  reducing  background  radiation:  spectral 
wave  filtering  technology.  This  method  reduces  background  photon  radiation,  radically  lowers 
the  system’s  noise  —  especially  photon  shot  noise  —  and  raises  the  system’s  signal-to-noise 
ratio. 

1  Artificial  satellite  spectrum  characteristics  and  background  characteristics 

1.1  Artificial  satellite  spectrum  characteristics  Artificial  satellites  are  not  usually  self- 
luminous.  The  sun  shines  on  them  in  the  upper  atmosphere  and  they  reflect  the  sun’s  light.  The 
intensity  of  illumination  of  spheroids  [viewed]  from  the  earth’s  surface  can  be  calculated  using 
the  following  formula: 


£(A)  =  r(A)seccz>  •  fi(A)  S£0(A) 

•  [(2  •  RV3  •  (Or  -  fi)  •  COS(/3)  +  SIN(fi))J/(r  •  L*)  (1) 
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Here,  Z  is  the  zenith  angle,  p(A.)  is  the  spectral  reflectivity  of  the  spheroid’s  surface,  R  is  the 
radius  of  the  spheroid,  (3  is  the  phase  angle  (the  included  angle  between  the  sun,  the  target,  and 
the  detector),  L  is  the  target’s  apparent  distance,  and  E0  is  the  sun’s  intensity  of  illumination 
outside  the  atmosphere. 
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The  surface  paint  colors  of  targets  being  detected  are  either  white  or  green.  The  spectral 
illumination  characteristics  (obtained  through  numerical  calculation)  of  targets  viewed  from  the 
earth  are  shown  in  Figure  1. 

1.2  Background  spectral  characteristics  During  the  day,  the  sky  is  bright.  Its 
luminance  ranges  from  0. 1  to  1  stilb111.  Its  spectral  irradiance  characteristics  are  shown  in  Figure 

1.  One  can  tell  from  Figure  1  that  background  radiation  is  primarily  concentrated  in  the  0.2  - 
0.6  pm  wave  band,  and  artificial  satellite  radiation  is  in  the  small  wave  band  from  0.2  -  0.4  pm. 
This  provides  conditions  for  the  application  of  spectral  wave  filtering  technology. 

2.  Analysis  of  spectral  filtering  technology  and  contrast 

Spectral  filtering  technology  is  the  choosing  of  an  appropriate  wave  filter  based  on  the 
difference  between  the  spectral  radiation  from  the  artificial  satellite  being  detected  and  the 
spectral  radiation  from  its  background.  It  causes  light  radiation  entering  the  photoelectric 
receiver  to  be  limited  to  necessary  wavelengths,  and  reduces  background  radiation  photons 
received  by  the  photoelectric  receiver  to  the  greatest  extent  possible,  thus  raising  target-signal 
contrast  and  signal-to-noise  ratio. 


The  definition  of  the  target  artificial  satellite’s  optical  contrast  Ct2]  is  the  ratio  between 
the  spectral  radiant  luminance  Ls  of  the  target  on  the  detector’s  surface  to  the  spectral  radiant 
luminance  of  the  background  Lg,  i.e.: 

C  =  L,/La  (2) 

It  can  be  proved  that  the  optical  contrast  given  in  the  above  definition  is  equivalent  to  that  in  the 
formulas  below: 


C  =  HJHa 


(3) 


r* 

Jit 

f 

.mi 


<A)  •  r*CA)  •  r,(X)dX 


E,tX)  •  r*(A)  •  rf(X)dX 


(4) 
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In  these  formulas,  Hs  and  HB  are,  respectively,  the  integral  illumination  of  the  artificial  satellite 
and  of  its  background;  r»  is  the  optical  system’s  spectral  penetration  rate;  Eg  and  Eg  are, 
respectively,  the  spectral  illumination  of  the  target  and  its  background  at  the  optical  system’s 
entrance  pupil;  and  tf  is  the  light  filter’s  spectral  penetration  rate. 

If  Ag  is  used  to  express  the  cutoff  wavelength  of  the  spectral  light  filter,  and  when  A  > 
Afl,  then  xF=  1,  and  when  A  <  Ag,  then  xF  =  0,  clearly,  choosing  a  light  filter  with  an  appropriate 
spectral  penetration  rate  according  to  differences  in  the  spectral  distribution  of  the  target  and  its 
background  causes  xF  to  have  different  effects  on  the  target  and  its  background.  Choosing  an 
appropriate  photoelectric  detector  allows  the  spectral  response  of  the  detector  to  match  the  light 
filter’s  penetration  rate  characteristics,  so  that  the  highest  amount  of  background  suppression  may 
be  attained.  In  this  way,  elevation  of  the  target’s  signal  contrast  is  achieved. 
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3,  The  relationship  between  contrast  and  signal-to-noise  ratio 

The  photoelectric  detectors  that  nonluminous  artificial  satellite  photoelectric  detection 
systems  can  use  have  CCD  area  array  devices  and  silicon  target-intensification  camera  tubes. 
Today,  most  use  CCD  area  array  devices.  Below,  the  detection  performance  of  photoelectric 
systems  is  calculated  using  CCD  area  array  detectors  as  photoelectric  receivers. 

3.1  Calculation  of  signal-to-noise  ratio  The  detection  performance  of  photoelectric 
systems  is  often  defined  as  the  signal-to-noise  ratio  of  the  system’s  output  signals.  Signal-to- 
noise  ratio  is  defined  as: 

SNR  =  Is/TK)l/i  (5) 

In  this  formula,  7S  represents  the  signal  current  of  a  CCD  pixel,  and  represents  the  mean- 
root-square  value  of  the  noise  current  of  a  CCD  pixel.  7S  can  be  written  as[41: 

Is  =  Ns  •  fH  .  e  (g) 

In  this  formula,  fH  represents  the  pixel’s  scanning  frequency,  e  is  the  electron’s  quantity  of 
electric  charge,  and  Ns  is  the  electric  charge  quantity  of  a  signal  output  by  a  CCD  pixel.  Ns  can 
be  obtained  using  the  following  formula141: 

Ns  —  Hs  'Qs'S'A'Ts*  K/e  (7) 
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In  this  formula,  A  is  the  system’s  entrance  pupil  point,  S  is  the  sensitivity  of  the  CCD,  Hs  is  the 
signal  light’s  integral  illumination,  Qs  is  the  quantum  efficiency  of  the  CCD  to  the  signal  light, 
Ts  is  the  signal  light’s  integral  time,  and  K  is  the  penetration  rate  of  a  neutral  light  filter. 

The  CCD’s  noise  sources  can  be  summed  up  as  follows:  shot  noise,  transfer  noise,  and 
thermal  noise. 

Thus,  (In)1/:  can  be  written  as: 

(  TZ  ',1/2  —  77"?  ’  TT  "TTy/I 

\*n)  —  UWJ  T  iNj  .  *Ni)  £g^ 

1.  Shot  noise 


The  randomness  of  incident  photon  current  produces  incident  photon  shot  noise.  The 
mean  square  current  of  shot  noise  produced  by  signal  current  7S,  background  current  /„,  and  dark 
current  /D  is  a  function  of  the  CCD  detector  band  A f.  Its  expression  is'5*: 


I2n\  —  2  •  e  •  (7,  -f-  It  +  J0)  .  A  f  (9) 

Because  each  pixel  of  the  CCD  can  be  seen  as  an  integral  wave  filter  with  integral  time 
bandwidth  istsl: 

A /  =  1/(2  •  T)  (10) 

The  background  current  from  a  pixel  of  a  CCD  can  be  written  as[4]: 

-  N,  *  fH  -e  <H) 

In  this  formula,  NB  is  the  background  quantity  of  electric  charge  from  a  CCD  pixel.  NB  c 
written  as[4]: 

Nb  =  Hb  •  Qb  •  S  •  A  •  Tb  •  Af  •  K/e  (22) 


In  this  formula,  HB  is  the  integral  intensity  of  background  light,  QB  is  the  quantum  efficiency  of 

the  CCD  to  background  light,  Ts  is  the  integral  time  of  background  light,  and  AP  is  the  angle  area 
of  a  CCD  pixel. 


Because  dark  current  /D  « IB,  when  we  do  not  consider  the  shot  noise  produced  by  dark 
current,  can  be  expressed  as: 


/»!=«•  (Ij  +  Ib)/T 


(13) 
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2.  Transfer  noise 


The  kind  of  transfer  noise  that  has  the  greatest  effect  is  transfer  loss  noise,  which  is  also 
a  kind  of  shot  noise.  It  is  primarily  produced  by  7S  and  /B.  The  mean  square  current  of  this  kind 
of  shot  noise  is[4]: 
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7 &-(«•/.  +  «•/.)•  2  tN/T  (14> 

where  eN  represents  the  overall  transfer  loss  rate  following  N  transfers.  In  general,  eN  =  0.1. 


3.  Thermal  noise 

Because  thermal  noise  is  much  fainter  than  shot  noise,  we  will  omit  thermal  noise  here. 


The  system’s  signal-to-noise  ratio  can  be  expressed  as: 


SNR  =  U(PN)m  = 


_ Hs  •  A  •  S  •  Ts  •  K  '  _ 

0  ♦  5  •  A  •  K  •  1.  2  (Hs  -Qs  -Ts  +  H,  -Q,  ‘Tb  - 


When  detecting  artificial  satellites  against  a  clear  daylight  background,  because  Hs«  HB  and 
Ts  <  Tb,  signal-to-noise  ratio  can  be  simplified  to: 


SNR  = 


Hs  •  Qs  •  Ts  -  (A  ♦  5  -  K}  1/2 
(.e  •  1.  2  •  •  Qt  •  T B  •  Af )  1/2 


(15) 


3.2  The  relationship  between  contrast  and  signal-to-noise  ratio  Assuming  that  the 
brightness  of  the  sky  background  is  definite,  from  formula  (3)  we  can  obtain: 

Hs=C-Ha 


and  because  TB  =  T,  therefore: 


SNR  = 


C  •  •  T.j  4  (.H b  *  S  *  A  m  1/2 

~e  ♦  1.  2  •  Qa  •  T  •  Ar)  1/2 


(17). 


When  the  background  brightness  detected  by  the  detector  surface  is  definite,  the  artificial 
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satellite’s  signal-to-noise  ratio  forms  a  direct  ratio  with  its  optical  contrast. 

4.  An  analysis  of  the  detection  capabilities  of  photoelectric  detection  systems 

Detection  capability  is  an  important  problem  that  must  be  solved  first  in  photoelectric 
detection  systems.  The  detection  capability  of  photoelectric  detection  systems  is  decided  by  the 
following  factors161: 

(1)  The  artificial  satellite’s  overall  dimension,  its  surface  reflectivity,  and  the  altitude  of 
its  orbit; 

(2)  The  zenith  angle,  phase  angle,  and  background  conditions  when  observing  the  artificial 
satellite; 

(3)  The  optical  system’s  parameters  (caliber,  penetration  rate,  relative  aperture,  etc.)  and 
the  spectral  filter’s  cutoff  wavelength; 

(4)  Functions  of  the  receiver  (spectrum  corresponds  to  quantum  efficiency,  saturation 
electron  number,  etc.). 

According  to  the  Poisson  distribution,  when  assured  that  SNR  =  6  in  95%  of  detection 
probability,  by  using  formula  (15)  we  can  obtain: 

„  _  6  •  (1.  2  •  *  •  H,  •  Q„  •  T  •  A,)  1/2 
s  Qs  •  T,(5  .  A  •  20 1/2  U8) 

After  finding  the  target’s  Hs,  it  can  be  changed  into  the  target’s  visual  magnitude  Mv.  The 
results  of  numerical  calculations  show  that  spectral  wave  filter  technology  is  effective  at  raising 
the  system’s  limits  of  detection.  Without  using  spectral  wave  technology,  where  SNR  =  6,  the 
CCD’s  saturation  electron  number  is  130  K,  telescope  diameter  D  =  400  mm,  background 
luminance  is  5.1  stellar  magnitudes,  and  the  artificial  satellite  surface  is  painted  white,  gray,  or 
black,  the  dimmest  detectable  artificial  satellite  is  6.6  stellar  magnitudes;  however,  when  a 
spectral  wave  filter  with  cutoff  wavelength  of  0.5  pm  is  added,  the  dimmest  detectable  artificial 
satellite  is  7.0  magnitudes.  When  the  artificial  satellite  surface  is  painted  green,  without  spectral 
wave  filter  technology,  the  dimmest  detectable  target  is  6.4  stellar  magnitudes;  but  when  a 
spectral  wave  filter  with  cutoff  wavelength  of  0.5  pm  is  added,  the  dimmest  detectable  target  is 
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6.9  stellar  magnitudes,  and  when  a  spectral  wave  filter  with  cutoff  wavelength  of  0.7  pm  is 
added,  the  dimmest  target  the  system  can  detect  is  7.1  stellar  magnitudes.  Clearly,  spectral  wave 
filter  technology  can  effectively  suppress  background  photon  radiation,  and  thus  raise  the 
system’s  detection  limits. 
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5.  Numerical  calculations 

Four  surface  paint  colors  were  used  as  examples:  white,  green,  gray,  and  blackm. 
Numerical  calculations  were  carried  out  based  on  their  surface  reflectivity  data.  The  results  are 
displayed  in  figures  2,  3,  and  4.  In  these  figures,  ^  represents  the  cutoff  wavelength  of  the 
spectral  wave  filter.  Conditions  where  =  0.2  pm  are  equivalent  to  conditions  where  spectral 
wave  filters  are  not  used.  We  can  tell  from  the  figures  that  no  matter  what  color  the  target  is 
painted,  using  spectral  wave  filter  technology  can  markedly  improve  the  contrast  and  signal-to- 
noise  ratio  of  the  target’s  signals,  and  can  raise  the  detection  limits  of  photoelectric  detection 
systems.  Calculation  results  are  shown  in  tables  1,  2,  and  3. 

6.  Conclusions 


Results  of  numerical  calculation  show  that  spectral  wave  filter  technology  is  effective  at 
enhancing  the  artificial-satellite  detection  capabilities  of  visible  light  detection  systems  in 
daylight.  When  detecting  nonluminous  artificial  satellites,  satellites’  optical  contrast  can  be 
multiplied  by  4  to  8  times,  and  their  signal- to-noise  ratio  can  be  multiplied  by  1.5  to  2.5  times. 
In  conditions  where  the  system’s  parameters  are  definite,  the  system’s  limit  detection  magnitude 
can  be  raised  by  0.4  to  0.7  visual  magnitude. 
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Key:  (1).  Table  1.  (2).  Wavelength  unit:  ym.  (3).  Contrast.  (4).  Painted  white.  (5).  Painted 
green.  (6).  Painted  gray.  (7).  Painted  black. 


Key:  (1).  Table  2.  (2).  Wavelength  unit:  ym.  (3).  Signal-to-noise  ratio.  (4).  Painted  white.  (5). 
Painted  green.  (6).  Painted  gray.  (7).  Painted  black. 
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Key:  (1).  Table  3.  (2).  Unit:  Visual  magnitude.  (3).  Detection  limit.  (4).  Painted  white.  (5). 
Painted  green.  (6).  Painted  gray.  (7).  Painted  black. 
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Key:  (1).  Figure  1.  Spectral  intensity  of  the  target  and  its  background.  (2).  Spectral  intensity. 
(3).  Background  characteristics.  (4).  Characteristics  of  the  target  when  painted  white.  (5). 
Characteristics  of  the  target  when  painted  green.  (6).  Figure  2.  The  relationship  between  contrast 
and  Xq.  (7).  Contrast.  (8).  Painted  green.  (9).  Painted  white.  (10).  Painted  gray.  (11).  Painted 
black. 
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Key:  (1).  Figure  3.  The  relationship  between  signal-to-noise  ratio  and  Aq.  (2).  Signal-to-noise 
ratio.  (3).  Painted  gray.  (4).  Painted  white.  (5).  Painted  green.  (6).  Painted  black.  (7).  Figure 
4.  The  relationship  between  detection  limits  and  V  (8).  Detection  limits.  (9).  Painted  gray. 
(10).  Painted  white.  (11).  Painted  green.  (12).  Painted  black. 
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